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An explicit expression is given for the correlation function of blinking systems, i.e. systems ex- 
hibiting light and dark periods in their fluorescence. We show through the example of terrylene in a 
crystalline host that it is possible to determine by means of this explicit expression photo-physical 
parameters, like Einstein coefficients and the mean light and dark periods by a simple fit. In addition 
we obtain further parameters like the frequency of the various intensity periods and the probability 
density of photons scattered off the host crystal. It turns out that this approach is simpler and 
allows greater accuracy than previous procedures. 
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I. INTRODUCTION 

Since about twenty years it has been possible to ob- 
serve and study the fluorescence of single ions in Paul 
traps and, more recently, also of single molecules em- 
bedded in a crystal In the presence of metastable 
triplet levels the molecule can exhibit light and dark pe- 
riods in its fluorescence ("blinking") la 0, just as ions 
in the or A configuration i7|,ia,liliaEllllllllIl- 
Spectroscopy of single molecules is of fundamental im- 
portance both for basic quantum mechanical aspects like 
antibunc hing, q uantum jumps or the dynamical Stark 
effect 0, EsllTal. as well as for applications in chemical 
and biophysics, cf. e.g. For these apphcations a 
precise knowledge is required of the parameters charac- 
teristic of fluorescence, like Einstein coefficients and the 
mean duration of intensity periods. Often these param- 
eters cannot be calculated directly but have to be deter- 
mined indirectly from suitable experimentally accessible 
quantities. 

One of the most important statistical quantities for the 
description of fluorescing quantum systems is the inten- 
sity correlation function, ^(t) Its behavior for small 
times T yields insight into the nature, classical or quan- 
tum, of the photon statistics i.e. bunching {g{0) > 1) or 
antibunching {g{0) < 1), respectively, and it shows the 
effect of Rabi oscillations. In addition it may yield in- 
formation on light and dark periods through long-time 
correlations. The importance of ^(t) originates from its 
easy experimental accessibility as it does not depend on 
the detector efficiency. To obtain a deeper physical un- 
derstanding of this important quantity, analytical results 
are often useful. By means of our recently proposed ap- 
proach for an analytical calculation of photon correlation 
functions of arbitrary blinking quantum systems jl9j we 
derive in this paper, for the first time, an explicit expres- 
sion for the correlation function of the four-level scheme 
used to describe the fluorescence of single molecular sys- 
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terns like terrylene or pentacene. In contrast to previous 
work, our result depends directly on the relevant photo- 
physical parameters, such as Einstein coefhcients, Rabi 
frequencies or the mean duration of the light and dark 
periods, with no further unknown quantities. This al- 
lows an easy determination of these in general unknown 
molecular parameters by a single fit to experimental data, 
where in comparison to previous approaches the fit is 
greatly facilitated through the explicit knowledge of ^(t), 
as will be shown in the present paper for the example of 
a terrylene molecule embedded in p-terphenyl. This in- 
tensity correlation function has recently been measured 
in Refs. [23, 0, 113 ■ In this paper we first derive an 
explicit expression of gir) for terrylene, taking into ac- 
count both its metastable levels, and then fit g{T) to the 
data of Ref. ^23] . From this we obtain the electronic Ein- 
stein coefficients, the mean light and dark periods and 
further fluorescence parameters and compare our results 
with those in the literature. 



II. THE INTENSITY CORRELATION 
FUNCTION 

The energy-level configuration of planar hydrocarbons, 
such as terrylene. can be described by an electronic three- 
level system [2J,|2J, with a singlet ground state, an 
excited singlet state, |3), and a triplet state, [2). Due 
to the interaction between the magnetic moments of the 
electron spins the metastable triplet state of terrylene 
splits into two metastable sub-levels with different pop- 
ulation and depopulation rates (zero-field splitting) j^. 
Each electronic level is characterized by a host of vibra- 
tional degrees of freedom (cf. Fig. 1(a)). These vibra- 
tional states have a life time of about 10~^^ s 1^. On the 
other hand, a transition from |3) to the vibrational states 
of |1) occur on a time scale of 10^^ s and from |3} to the 
triplet state |2) about every 10"'^ s. Therefore these tran- 
sitions can be described by effective electronic Einstein 

coefficients ^131,^32 ,^2*1 = li2, and by the Rabi fre- 
quency 1^31 of the laser driving the |1) ^ |3) transition, 
while the vibrational levels can be neglected. Hence the 
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dark periods to the light period by ^^id ■ These 
transition rates can be calculated in a fully quantum me- 
chanical way by the methods of Ref. 
given by (cf. Appendix A) 
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A fluorescence trajectory of terrylene is a Markovian 
jump process with three intensity steps, i.e. light, dark 
D^-*^^ and dark D*^^). Recently, the present authors 
have derived a highly accurate expression for the inten- 
sity correlation function of arbitrary blinking systems 
which is given in Appendix B and takes here the form 
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FIG. 1: Level scheme of terrylene. (a) Singlet ground 
state excited singlet state |3), and triplet state |2), as- 
sociated with a host of vibrational states. ISC (intersys- 
tem crossing, singlet-triplet transition) with low probability. 
Zero-field splitting of the metastable triplet state into two 
sublevels with different population and depopulation rates, 
(b) The extremely fast relaxation of the vibrational levels 
leads to an effective four-level system with two metastable 
states, 12^'^-') and 12^^'), and effective Einstein coefficients 
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level configuration effectively consists of a four-level sys- 
tem as in Fig. 1(b), with two metastable states |2'^^)) and 
12^^^), and one has 
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A fluorescence trajectory of such a system consists of 
light and dark periods. During a light period, the sub- 
system {|1),|3)} behaves like a two-level system whose 
intensity is given by p5| 
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A transition from |3) to one of the metastable states |2'^^^) 
or 12^^^) each initiates a dark period (denoted by D^^^ and 

D^^^) of mean duration A'^-^^ and A'^^i , respectively. The 
transition rates from the light period to the two dark pe- 
riods are denoted by and Plq, and those from the 



5(t) = Pll{t) g2ir) , 
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where Pll (t) is the conditional probability to have a light 
period at time t = t provided that at i = there was 
also a light period, where 



Pl = lim Pll(t) 



(6) 



is the probability to flnd a light period at all, and where 
g2{T) is the intensity correlation function of the {|1), |3)} 
two-level subsystem. From Ref. one has 
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To determine Pll(t) we consider the conditional prob- 
ability, Plq(t), that at time t = t one has a dark pe- 
riod D*^*'' under the condition of a light period at i = 
{i = 1,2). These conditional probabilities satisfy the rate 
equations 
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which is easily seen by noting that the first terms on the 
r.h.s. describe the decrease of the probabilities and the 
second their increase. This can easily be generalized for 
an arbitrary number of light and dark periods as shown 
in Appendix B. With the initial condition Pll(O) = 1 
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one finds in a straightforward way 
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From this and from Eq. lO one obtains 
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Thus the correlation fmiction ^(t) in Eq. Q f^or the four- 
level system of Fig. 1(b) satisfying the inequalities in 
Eq. is given in terms of the photo-physical param- 

eters A31 , f^gi , ^^2^ , a'§ , ^ , and A^^^ . 



One can re-express gir) through the mean durations 
of the three periods, denoted by TliTq^-* and Tj^\ re- 
spectively, and by the branching ratios 
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With these relations one obtains, after some calculation, 
the explicit expression 
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where 
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r = 2 - + (1 - 2p,)/nr + Ap,p,/T, 

(18) 

One sees that the statistics of intensity periods leads to 
the well known bi-exponential tail in the intensity corre- 
lation function at a much larger time scale than where 
g2iT) plays a role. Therefore one has a hump at interme- 
diate times r ~ 10~^s, where g2{T) ~ 1 and Pll(t) ~ 1. 
The height of this hump is given by 1/Pl, according to 
Eq. © and to Eq. ^ for r ^ and 32(1") ~ 1. More- 
over, with the explicit form of ^(r) in Eq. (fT7|l . all photo- 
physical parameters can be obtained by a single fit to 
experimental data, as shown in the following. 



III. EXAMPLE: FIT AND RESULTS FOR 
TERRYLENE 



We use the experimental data for the intensity correla- 
tion function of terrylene obtained in Ref. 22J . Details of 
the experimental setup can be found in Ref. Since 
some of the laser light is scattered off the embedding crys- 
tal, with intensity I^c , say, this gives rise to additional 
poissonian correlations. Therefore the intensity correla- 
tion function in a light period is modified to the weighted 
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with /l given by Eq. jSJl. This has to be inserted into 
Eq. l(T7|) instead of 52 (t). 

The thus modified Eq. H17() has been fitted to the ex- 
perimental data. In Eq. H17|l there are two different time 
scales - for times r > 10~^s the first factor, 5™°'^(t) is 
essentially 1 and for times r < 10~^ s the second factor 
is practically constant. Therefore, we have first fitted 
the second factor for r > 10~^s and then fitted .g™°'*(T) 
multiplied by a constant for r < 10^^ s. Thus, in each of 
the two fits one has at most four free parameters, which 
reduces the numerical error. The fitted curve and the 
experimental data are shown in Fig. 2, and the resulting 
parameter values are listed in Table I. 

With the results of Table I, one now could calculate 



and A^2i from Eqs. ©, H, and d - (O. How- 
ever, numerically it is better to insert the expressions for 
ThjTj^^ and pi into Eq. (|17|l and then again fit the sec- 
ond factor to the experimental data for r > 10^^ s. With 
this procedure one obtains all Einstein coefficients of the 
intersystem crossings, as listed in Table IL 
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TABLE I: Photo-physical parameters for terrylene in p-terphenyl, obtained by a fit of Eq. 1171 to the data of Ref. and 
compared with resuhs in the literature. 
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Ref. [22j 
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FIG. 2: The intensity correlation function of a single fluo- 
rescent terrylene molecule embedded in a host crystal of p- 
terphenyl. Shown are the experimental data from Ref. |23l 
(x) and the fit function obtained with Eq. p7^ (solid line). 
This function contains both the Rabi oscillations for small 
times as well as the long time correlations (hump) due to the 
metastable levels. 



Our value for A31 agrees with that of Refs. [IS US] within 
the error bars. 

In case of the Einstein coefficients of the intersystem 
crossin gs o ur values quahtatively agree with those of 
Ref. Hillll, where the results of Refs. depend 
on the same data as ours. There is overall agreement in so 
far as the population rates of the metastable states 
are an order of magnitude smaller than the corresponding 
depopulation rates and that the values for i = 1 are 
an order of magnitude smaller than those for i — 2. As a 
consequence, the metastable state |2(^)) is occupied con- 
siderably less often than 12^^^}, a fact also mirrored in the 
largely different branching ratios pi and p2- Hence the 
dark period associated with the state |2^^)} occurs much 
less frequently than the other one, and hence a statisti- 
cal determination may be more error prone than our fit 
procedure. In Table II, somewhat larger differences are 
seen for A''2i , the population rate of the less frequently 
occupied metastable state. 



IV. CONCLUSIONS AND OUTLOOK 



TABLE 11: Einstein coefficients obtained here for the singlet- 
triplet transitions (ISC), compared with values in the litera- 
ture. 
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Ref. [27] 
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Ref. [22j 
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Ref. [5] 
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The mean durations of the light and dark periods ob- 
tained here essentially agree with those of Ref. 22, 27], 
with a difference of about 23% in the case of Tl and 
with a drastically reduced error bar in our value for Tj^\ 
The reason probably is that a direct determination of the 
period duration as done in Refs. |23,|23| depends sensi- 
tively on the choice of the time window employed [26j |. 



In this paper an explicit expression for the intensity 
correlation function of blinking systems has been applied 
to a four-level system with two metastable states. Us- 
ing terrylene in p-terphenyl as an example it has been 
shown that the relevant photo-physical parameters can 
be obtained from a single fit to experimental data for 
the correlation function. This is a considerable simplifi- 
cation compared to previous procedures, and it can also 
result in smaller error bars. In this way we have ob- 
tained the Einstein coefficient of the fast transition, the 
mean durations of the light and dark periods, as well as 
the population and depopulation rates of the metastable 
levels. Our results qualitatively agree with values given 
in the literature. In addition we have obtained further 
parameters like the Rabi frequency fl^i and the contri- 
bution, /sc, to the correlation function due to the light 
back-scattered off the crystal. 

It should also be possible to apply our approach of 
an explicit calculation of the correlation function not 
only to terrylene but also to other blinking systems with 
a more complicated statistics of light and dark peri- 



5 



ods. Particularly interesting in this context would be 
the recently performed fluorescence measurement of two 
dipole-interacting molecules and to determine the 
interaction parameters of such a system by the present 
methods. 



APPENDIX A 

We indicate here the derivation of Eqs. Q and Q 
for the transition rates between the different fluorescence 
periods of the level system in Fig. 1(b). To do this we 
use the method of Ref. [2^. During a dark period D^*^ 
the system is, to high accuracy, in the state |2(*)), with 
corresponding density matrix 



p« = |2«)(2W| 



(Al) 



and during a light period in the steady state of the 
{|1), |3)} subsystem, with the density matrix 



3j;^((4. + o^,)|i)(il 
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+nl,\3){3\+iA3,n3,{\l){3\ - |3)(1|) . (A2) 
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Starting from one of these density matrices, the increase 
in occupation of the other subsystems during a time At 
yields the transition rates; here At has to satisfy 

(42)-', (At)-' « At « A^,\n^,' . (A3) 
With initial state p{to) = pi^ one then has 
d 



p« =-(2W|p(t)|2W),=,„+A* 



and for p(to) — p^^ one has 



(A4) 



PdL = :t7 (3|p|3) + (l|p|l) ^ . (A5) 
at\ / t=to+At 

To evaluate this we use the Bloch equations which can 
be written in the form j2^ ^] 

^p(t) ^ ^^(^Hcp{t)-p{t)Hl)+n{p{t)) (A6) 
= Lp{t) , (A7) 
where the so-called conditional Hamiltonians He is given 

by 

He = |(^|3)(3|+4V|2«)(2W|+4i'|2(^)>(2(^)|) 

+ ^(|1)(3| + |3)(1|) (A8) 
and the super-operator TZ by 

2 

TZ{p) = 5]|l)(^3i(3|p|3)+4?(2«|p|2«))(l| 

1=1 

+4!)|2«)(3|p|3)(2«| . (A9) 



Using Eq. IjAGp one can show that Eqs. HA4|) and (|A5|I 

can be written as 

P'li - 42(3|p(to + At)|3)-A«(2W|p(to + At)|2«) 

(AlO) 

p'dI - E(^2i^(2(")|p(to-l-At)|2(")) 



a=l 

-4^)(3|p(to + At)|3) 



(All) 



The necessary matrix elements of p{to + At) are obtained 
from the Bloch equations written as in of Eq. HA7|I . Writ- 
ing L in the obvious form 



L = Lo(A31,f^3l)+Li(42-4?) 



(A12) 



one obtains in a way analogous to the usual quantum 
mechanical perturbation theory, to first order in Ag!] and 



p{to + At) 



e^'^'pito) 
e^°^V(to) 



/ dr e^°('^*~")Lie^«Xto) . 
Jo 

(A13) 



Since p(to) is a stationary state of subsystems one has 
Lop(to) = and hence e^°^ p(to) = p{to). After a change 
of variable one thus obtains 



At 



p{to + At) = p{to) + I dr e^'>^Up{to) . (A14) 

JQ 

Now Lip(to) can be decomposed into a component par- 
allel to the null space of Lq and a component belonging 
to the subspace for nonzero eigenvalues, 



Lip(io) = P||Lip(io) + (1 - P|| )Lip(to) . (A15) 



Here 



P|l = — i dz(z-Lo)-i , (A16) 
27ri /c, 

where Co is a path in the complex plane enclosing no 
other eigenvalue than zero ^l|. The integral can be cal- 
culated by the residue method. In view of Eq. ljA3|l the 
contribution to Eq. HA14|I of the P|| term in Eq. (|A15|) is 
neghgible. The contribution to Eq. (|A14p of the (1 -Py) 
term is governed by the nonzero eigenvalues of Lq which 
possess a (large) negative real part of the order of A^i 
and r^si . Therefore the integrand in Eq. IIA14|) is strongly 
damped in r, and in view of Eq. ljA3p the upper integra- 
tion limit can be replaced by oo, yielding 

p{to + At) - p{to)+ / dTeI^«^(]l-P[|)Lip(io) 



p(to) + (e-Lo)-i(l-P||)Lip(to) 



(A17) 
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with e +0. Inserting this into Eqs. (|Alip and (jA10|) 
one obtains, after some calculations involving 16x16 ma- 
trices, the transition rates of Eqs. (PJ and l@J. These hold 



to first order in Ai^} and A 



32 ■ 



APPENDIX B 

We give here generalized expressions for the intensity 
correlation function of fluorescing systems with an arbi- 
trary number of different intensity periods with intensity 
li . A detailed discussion can be found in Ref. [T9l |. 

Let Pi be the probability for the occurrence of period 
i and let Pij (r) be the probability to have period j at 
time r provided one had period i at r = 0. Then it 
can be shown, that a highly accurate expression for the 
normalized intensity correlation function is 



(Bl) 



where gi (r) is the correlation function within a given pe- 
riod i and usually easier to calculate than that of the 
complete system. In the case of only one light period 
and two dark periods Eq. (|Bip reduces to Eq. ((S)). 

Let Pij be the transition rates from period i to j. Then 
the Pijir) are easily seen to obey rate equations, e.g. 



Al(r) = [-J2pik)Pllir)+P2lPl2iT) + ---+P„lPlniT) . 

k 

(B2) 

In general, with the matrix B = {Bij), 

Bij = Pij - Sij ^Pik , (B3) 



and the matrix 



one has 



(B4) 



P = PB , (B5) 
with the initial condition Pij (0) — Sij , or 

P(0) = 1 . (B6) 



The solution of equation (jB5|l with this initial condition 
can be written as 



P(r) 



(B7) 



If /io, • • • , Mn-i are the eigenvalues of B (assumed dis- 
tinct) then _32J 



n-l 



n / - Mi Met 

The properties of the matrix B are closely related to 
those of stochastic matrices 's^ , and under quite general 
conditions B has a single eigenvalue = and eigen- 
values /ii, . . . , with negative real part. 

To find the Pi's, we note that for t — + oo the memory 
to the initial condition is in general lost. Therefore, for 
any k. 



(B9) 



[1] W. Neuhauser, M. Hohenstett, P. Toschek, and 

H. Dehmelt, Phys. Rev. A 22, 1137 (1980). 
[2] W.E. Moerner and L. Kador, Phys. Rev. Lett. 62, 2535 

(1989) . 

[3] M. Orrit and J. Bernard, Phys. Rev. Lett. 65, 2716 

(1990) . 

[4] P. Tamarat, A. Maali, B. Lounis, and M. Orrit, J. Phys. 

Chem. A 104, 1 (2000). 
[5] T. Basche, S. Kummer, and C. Brauchle, Nature 373, 

132 (1995). 

[6] M. Vogel, A. Gruber, J. Wrachtrup, and C. von Bor- 
czyskowski, J. Chem. Phys. 99, 14915 (1995). 

[7] H.G. Dehmelt, Bull. Am. Phys. Soc. 20, 60 (1975) 20, 
60 (1975). 

[8] W. Nagourney, J. Sandberg, and H. Dehmelt, Phys. Rev. 

Lett. 26, 2797 (1986). 
[9] J.C. Bergquist, R.G. Hulet, W.M. Itano, and 
D.J. Wineland, Phys. Rev. Lett. 57, 1699 (1986). 
[10] T. Sauter, W. Neuhauser, R. Blatt, and P.E. Toschek, 

Phys. Rev. Lett. 57, 1696 (1986). 
[11] T. Sauter, R. Blatt, W. Neuhauser, and P.E. Toschek, 
Opt. Comm. 60, 287 (1986). 



[12] F. Diedrich and H. Wahher, Phys. Rev. Lett. 58, 203 
(1987). 

[13] W.M. Itano, J.C. Bergquist, R.G. Hulet, and 

D.J. Wineland, Phys. Rev. Lett. 59, 2732 (1987). 
[14] W.M. Itano, J.C. Bergquist, and D.J. Wineland, Phys. 

Rev. A 38, 559 (1988). 
[15] T. Basche, W.E. Moerner, M. Orrit, and H. Talon, Phys. 

Rev. Lett. 69, 1516 (1992). 
[16] P. Tamarat, B. Lounis, J. Bernard, and M. Orrit, Phys. 

Rev. Lett. 75, 1514 (1995). 
[17] P. Schwille, Cell. Biochem. Biophys. 34, 383 (2001). 
[18] L. Mandel and E. Wolf, Optical coherence and quantum 

optics (Cambridge University Press, Cambridge, 1995). 
[19] G.C. Hegerfeldt and D. Seidel, Quantum Semiclass. Opt. 

4, 245 (2002). 

[20] S. Kummer, T. Basche, and C. Brauchle, Chem. Phys. 

Lett. 229, 309 (1994). 
[21] S. Kummer, S. Mais, and T. Basche, J. Phys. Chem. 99, 

17078 (1995). 

[22] T. Basche and C. Brauchle, Ber. Bunsenges. Phys. Chem. 

100, 1269 (1996). 
[23] J.B. Birks, Organic molecular photophysics (Wiley, 



7 



1973). 

[24] J. Bernard, L. Fleury, H. Talon, and M. Orrit, J. Chcm. 

Phys. 98, 850 (1993). 
[25] H. Carmichael and D.F. Walls, J. Phys. B 9, 1199 (1976). 
[26] S.U. Addicks, A. Beige, M. Dakna, and G.C. Hegerfeldt, 

Eur. Phys. J. D 15, 393 (2001). 
[27] S. Kummer, Ph.D. thesis, Ludwig-Maximilians- 

University Miinchon (1996). 
[28] C. Hettich, private communication. 



[29] G.C. Hegerfeldt, Phys. Rev. A 47, 449 (1993). 

[30] G.C. Hegerfeldt and D.G. Sonderrnami, Quantum Semi- 
class. Opt. 8, 121 (1996). 

[31] T. Kato, Perturbation theory for linear operators 
(Springer, Berlin, 1966). 

[32] F. Gantmacher, The theory of matrices (Chelsea Publ. 
Co., 1959). 



